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The gravastar picture is an alternative model to the concept of a black hole, where there is an 
effective phase transition at or near where the event horizon is expected to form, and the interior 
is replaced by a de Sitter condensate. In this work, a generalization of the gravastar picture is 
explored, by considering a matching of an interior solution governed by the dark energy equation of 
state, uj = p/p < — f/3, to an exterior Schwarzschild vacuum solution at a junction interface. The 
motivation for implementing this generalization arises from the fact that recent observations have 
confirmed an accelerated cosmic expansion, for which dark energy is a possible candidate. Several 
relativistic dark energy stellar configurations are analyzed by imposing specific choices for the mass 
function. The first case considered is that of a constant energy density, and the second choice, 
that of a monotonic decreasing energy density in the star's interior. The dynamical stability of the 
transition layer of these dark energy stars to linearized spherically symmetric radial perturbations 
about static equilibrium solutions is also explored. It is found that large stability regions exist 
that are sufficiently close to where the event horizon is expected to form, so that it would be 
difficult to distinguish the exterior geometry of the dark energy stars, analyzed in this work, from 
an astrophysical black hole. 
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I. INTRODUCTION 



The structure of relativistic stars and the phenomenon 
of gravitational collapse is of fundamental importance in 
astrophysics and has attracted much attention in the rel- 
ativist community since the formulation of general rela- 
tivity. Relatively to the construction of theoretical mod- 
els of relativistic stars, one may refer to the pioneering 
work of Schwarzschild , Tolman , and Oppenheimer 
and VolkofF '3| . Schwarzschild considered analytic so- 
lutions describing a star of uniform density, and Tolman 
developed a method providing explicit solutions for 
static spheres of fluid, which proved important to the 
study of stellar structure. Oppenheimer and VolkofF 3] 
by suitably choosing specific Tolman solutions studied 
the gravitational equilibrium of neutron stars, using the 
equation of state for a cold Fermi gas, consequently laying 
down the foundations of the general relativistic theory of 
stellar structures. Paging through history one finds an 
extremely extensive literature however, one may re- 
fer to the important contributions of Chandrasekhar 5] 
in the construction of white-dwarf models by taking into 
account special relativistic effects in the electron degener- 
acy equation of state, where he discovered that no white 
dwarf may have a mass greater than ~ 1.2 solar masses, 
which has been called the Chandrasekhar limit; and the 
work of Baade and Zwicky (6j, where they invented the 
concept of a neutron star, and identified astronomical ob- 
jects denoted as supernovae, representing a transitional 
collapse of an ordinary star into a neutron star. Far from 
undertaking an exhaustive review, one may also refer to 
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the work of Wyman, where in Ref. (2j isotropic coor- 
dinates were used to solve the relativistic equations of a 
perfect fluid with constant energy density, and in Ref. '3| , 
a critical examination and generalization of the Tolman 
solutions was undertaken; Buchdahl |^ and Bondi |l3 
generalized the interior Schwarzschild solution to more 
general static fluid spheres in the form of inequalities 
involving the mass concentration, central energy density 
and central pressure; Leibovitz also generalized some 
of Tolman's solutions by applying a more physical ap- 
proach; and the discovery of pulsars by Hewish et al |l2j| . 
the idea advanced by Gold 113] that pulsars might be 
rotating neutron stars, which was later confirmed by ob- 
servations, and the idealized models of rotating neutron 
stars 01 . It is also interesting to note that Bayin found 
new solutions for static fluid spheres 0|, explored the 
time-dependent field equations for radiating fluid spheres 
|T6| . and further generalized the analysis to anisotropic 
fluid spheres jlTj . 

Relatively to the issue of gravitational collapse, before 
the mid-1960s, the object now known as a black hole, 
was referred to as a collapsed star j!^]. Oppenheimer 
and Snyder 0|, in 1939, provided the first insights of 
the gravitational collapse into a black hole, however, it 
was only in 1965 that marked an era of intensive research 
into black hole physics. Although evidence for the ex- 
istence of black holes is very convincing, it has recently 
been argued that the observational data can provide very 
strong arguments in favor of the existence of event hori- 
zons but cannot fundamentally prove it j20t. This scep- 
ticism has inspired new and fascinating ideas. In this 
line of thought, it is interesting to note that a new fi- 
nal state of gravitational collapse has been proposed by 
Mazur and Mottola [2lj. In this model, and in the re- 
lated picture developed by Laughlin et al 0| , the quan- 
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turn vacuum undergoes a phase transition at or near the 
location where the event horizon is expected to form. 
The model denoted as a gravastar ((^raiitational vacuum 
star), consists of a compact object with an interior de 
Sitter condensate, governed by an equation of state given 
hy p — — p, matched to a shell of finite thickness with an 
equation of state p — p. The latter is then matched to an 
exterior Schwarzschild vacuum solution. The thick shell 
replaces both the de Sitter and the Schwarzschild hori- 
zons, therefore, this gravastar model has no singularity at 
the origin and no event horizon, as its rigid surface is lo- 
cated at a radius slightly greater than the Schwarzschild 
radius. It has been argued that there is no way of dis- 
tinguishing a Schwarzschild black hole from a gravastar 
from observational data '2^. It was also further shown 
by Mazur and Mottola that gravastars are thermodynam- 
ically stable. Related models, analyzed in a different con- 
text have also been considered by Dymnikova 23] . How- 
ever, in a simplified model of the Mazur-Mottola pic- 
ture, Visser and Wiltshire constructed a model by 
matching an interior solution with p = —p to an exterior 
Schwarzschild solution at a junction interface, comprising 
of a thin shell. The dynamic stability was then analyzed, 
and it was found that some physically reasonable sta- 
ble equations of state for the transition layer exist. In 
Ref. ^3 J ^ generalized class of similar gravastar mod- 
els that exhibit a continuous pressure profile, without 
the presence of thin shells was analyzed. It was found 
that the presence of anisotropic pressures are unavoid- 
able, and the TOV equation was used to place constraints 
on the anisotropic parameter. It was also found that the 
transverse pressures permit a higher compactness than 
is given by the Buchdahl-Bondi bound 9, 10] for perfect 
fluid stars, and several features of the anisotropic equa- 
tion of state were explored. In Ref. [Iflj, motivated by 
low energy string theory, an alternative model was con- 
structed by replacing the de Sitter regime with an interior 
solution governed by a Chaplygin gas equation of state, 
interpreted as a Born-Infeld phantom gravastar. 

It has also been recently proposed by Chapline that 
this new emerging picture consisting of a compact ob- 
ject resembling ordinary spacetime, in which the vacuum 
energy is much larger than the cosmological vacuum en- 
ergy, has been denoted as a "dark energy star" [23]. In 
fact, in the present paper a mathematical model general- 
izing the Mazur-Mottola picture, or for that matter, the 
Visser- Wiltshire model is proposed, where the interior de 
Sitter solution is replaced by a solution governed by the 
dark energy equation of state, p = ujp with oj < —1/3, 
matched to an exterior vacuum Schwarzschild solution. 
Note that the particular case oi oj — ~1 reduces to the 
Visser- Wiltshire gravastar model. The motivation for im- 
plementing this generalization comes from the fact that 
recent observations have confirmed that the Universe is 
undergoing a phase of accelerated expansion. Evidence 
of this cosmological expansion has been shown indepen- 
dently from measurements of supernovae of type la (SNe 
la) p8i] and from cosmic microwave background radiation 



|29| . A possible candidate proposed for this cosmic ac- 
celeration is precisely that of dark energy, a cosmic fluid 
parametrized by an equation of state uj = p/p < —1/3, 
where p is the spatially homogeneous pressure and p the 
dark energy density. If tu < —1, a case certainly not 
excluded, and in fact favored, by observations, the null 
energy condition is violated. For this case, the cosmic 
fluid is denoted phantom energy, and possesses peculiar 
properties, such as negative temperatures [s^] and the 
energy density increases to infinity in a finite time, re- 
sulting in a Big Rip [^J- It also provides one with a 
natural scenario for the existence of exotic geometries 
such as wormholes [sj Is^]- In this context, it was also 
shown that the masses of all black holes tend to zero as 
the phantom energy universe approaches the Big Rip |34| . 
However, it is interesting that recent fits to supernovae, 
CMB and weak gravitational lensing data favor an equa- 
tion of state with a dark energy parameter crossing the 
phantom divide lo = —1 ,35.. .36] . In a cosmological set- 
ting the transition into the phantom regime, for a single 
scalar field is probably physically implausible, thus 
the stress energy tensor should include a mixture of var- 
ious interacting non-ideal fiuids. 

As emphasized in Refs. [^S, IH^i in a rather different 
context, a subtlety needs to be pointed out: The notion 
of dark energy is that of a spatially homogeneous cos- 
mic fluid, however, it can be extended to inhomogeneous 
spherically symmetric spacetimes by regarding that the 
pressure in the dark energy equation of state is a neg- 
ative radial pressure, and the transverse pressure may 
be determined via the field equations. (An inhomoge- 
neous spherically symmetric dark energy scalar field was 
also considered in Ref. 37]). In this context, the gen- 
eralization of the gravastar picture with the inclusion of 
an interior solution governed by the equation of state, 
p — Lup with u; < —1/3, will be denoted by a dark energy 
gravastar, or simply a "dark energy star" in agreement 
with the Chapline definition We shall explore sev- 
eral configurations, by imposing specific choices for the 
mass function. We shall then explore the dynamical sta- 
bility of the transition layer of these models to linearized 
perturbations around static solutions, by applying the 
general stability formalism developed in Ref. 6S\, and 
which was also applied in the context of the stability of 
phantom wormholes js^ , and further analyze the evolu- 
tion identity to extract some physical insight regarding 
the pressure balance equation across the junction inter- 
face. The dark energy star outlined in this paper may 
possibly have an origin in a density fluctuation in the 
cosmological background. It is uncertain how such inho- 
mogeneities in the dark energy may be formed. However, 
a possible explanation may be inferred from Ref. [40| . 
where the dark energy equation of state was generalized 
to include an inhomogeneous Hubble parameter depen- 
dent term, possibly resulting in the nucleation of a dark 
energy star through a density perturbation. 

This paper is outlined in the following manner. In 
Section II, we present the structure equations of dark 
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energy stellar models. In Section III, specific models are 
then analyzed by imposing particular choices for the mass 
function. In Section IV, the linearized stability analysis 
procedure is briefly outlined, and the stability regions 
of the transition layer of specific dark energy stars are 
determined. Finally in Section V, we conclude. 



II. DARK ENERGY STARS: EQUATIONS OF 
STRUCTURE 



Consider the interior spacetime, without a loss of gen- 
erality, given by the following metric, in curvature coor- 
dinates 



-exp 

-r^ {dO 



g{f)df 
j'^ed(l)^), 



dr 



dr'^ 



1 — 2m(r)/r 



(1) 



where g{r) and m(r) are arbitrary functions of the radial 
coordinate, r. The function m(r) is the quasi-local mass, 
and is denoted as the mass function. The factor g{r) is 
the "gravity profile" and is related to the locally mea- 
sured acceleration due to gravity, through the following 
relationship: A = a/1 — 2m{r)/r g{r) |4ll. li^. The con- 
vention used is that g{r) is positive for an inwardly grav- 
itational attraction, and negative for an outward gravi- 
tational repulsion. Note that equivalently one may con- 
sider a function ^(r), defined as <I>(r) = — g{'r)dr, 
and denoted as the redshift function, as it is related to 
the gravitational redshift • 

The stress-energy tensor for an anisotropic distribution 
of matter is provided by 



T^,l, = ip + Pt)U^U^+ptgf,„ + [pr - PtjXi^Xu , 
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where C/^ is the four-velocity, is the unit spac elike 
vector in the radial direction, i.e., x^ = ^/l — 2m/r5^r- 
p{r) is the energy density, Pr{r) is the radial pressure 
measured in the direction of x'^j and pt(r) is the trans- 
verse pressure measured in the orthogonal direction to 
X^. 

Thus, the Einstein field equation, G^^, = SttT^^, where 
Gf^iy is the Einstein tensor, provides the following rela- 
tionships 
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where the prime denotes a derivative with respect to 
the radial coordinate, r. Equation |SJ| corresponds 
to the anisotropic pressure Tolman-Oppenheimer-Volkoff 
(TOV) equation. 

Now, using the dark energy equation of state, pr = cup, 
and taking into account Eqs. ||2Jl and we have the 



following relationship 
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Ther e is. however, a subtle point that needs to be empha- 
sized [Sj, 123 • The notion of dark energy is that of a spa- 
tially homogeneous cosmic fluid. Nevertheless, it can be 
extended to inhomogeneous spherically symmetric space- 
times, by regarding that the pressure in the equation of 
state p = ujp is a radial pressure, and that the transverse 
pressure may be obtained from Eq. JSJ. Note that for 
the particular case of w = —1, from Eq. 10, one has the 
following solution gtt — —(1 — 2m/r), which reduces to 
the specific class of solutions analyzed in Ref. . 
Using the dark energy equation of state Pr = ojp, Eq. 
in terms of the principal pressures, takes the form 

, fl + uj\ m + ujrm' 2, , , , 

Vr = -Vv —7 + - - Vr) , (7 

\ uj J r [r — 2m) r 

which taking into account Eq. (PJ, may be expressed in 
the following equivalent form 



A = 



87rr2 



2m' 



m rg 



(8) 



A = Pt — is denoted as the anisotropy factor, as it is a 
measure of the pressure anisotropy of the fluid comprising 
the dark energy star. A = corresponds to the particular 
case of an isotropic pressure dark energy star. Note that 
A/r represents a force due to the anisotropic nature of 
the stellar model, which is repulsive, i.e., being outward 
directed if pt > Pr, and attractive if pt < Pr- 

One now has at hand four equations, namely, the 
field equations ®-© and Eq. ©, with five unknown 
functions of r, i.e., p{r), Pr{r), Pt{r), g{r) and mlr). 
Obtaining explicit solutions to the Einstein fleld equa- 
tions is extremely difficult due to the nonlinearity of the 
equations, although the problem is mathematically well- 
defined. However, in the spirit of Ref. .33i| . we shall adopt 
the approach in which a specific choice for a physically 
reasonable mass function m{r) is provided and through 
Eq. (0, gir) is determined, thus consequently providing 
explicit expressions for the stress-energy tensor compo- 
nents. In the specific cases that follow, we shall consider 
that the energy density be positive and finite at all points 
in the interior of the dark energy star. 



III. SPECIFIC MODELS 

A. Constant energy density 

Consider the specific case of a constant energy density, 
p[r) = po, so that Eq. ^ provides the following mass 
function 



f ^ 47rpo 3 
m[r) — — - — r . 



(9) 
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Thus, using Eq. JBJ, one finds that g(r) is given by 

Ar{l + 3lu) 



1 - 2^r2 



(10) 



where for simphcity, the definition A — Anpo/S is used. 
Note that for uj < —1/3, we have an outward gravita- 
tional repulsion, g{r) < 0, which is to be expected in 
gravastar models. 

The spacetime metric for this solution takes the fol- 
lowing form 



1 - 2Ar^ 



[de^ 



sm 



(11) 



The stress-energy tensor components are given by Pr 
Lupo and 



Pt = wpo 



{1 + oj){l + 3uj)Ar'^ 
2w(l - 2Ar2) 



(12) 



The anisotropy factor is provided by 



Q 42 2 

A(i+^)(l + 3^). 



l-2Ar 



(13) 



One readily verifies that A<Ofor— l<a'<— 1/3 and 
A > in the phantom regime, cj < — 1. However, it is 
perhaps instructive to plot A, which is depicted in Fig. 
n Note that A = at the origin, r = 0, as was to be 
expected. For uj — ~1, then A = is also verified for 
arbitrary r. This latter condition is also readily derived 
from Eq. |(SJ), where taking into account pr — top for 
Pr — const, one verifies pt = Pr- 




FIG. 1: Plot of the anisotropy factor for a dark energy star 
with a constant energy density. We have defined the following 
dimensionless parameters: 5 — IS./ A and a — \J~Ar. We verify 
that A < in the range —1/3 > uj > —1 and A > in the 
phantom regime, ui < —1. 



B. Tolman-Matese- Whitman mass function 

Consider the following choice for the mass function, 
given by 

™(^)-2(rw)' (14) 

where bg is a non-negative constant. The latter may be 
determined from the regularity conditions and the finite 
character of the energy density at the origin r = 0, and 
is given by bg = Snpc/S, where pc is the energy density 
at r = 0. 

This choice of the mass function represents a mono- 
tonic decreasing energy density in the star interior, and 
was used previously in the analysis of isotropic fluid 
spheres by Matese and Whitman 44^ as a specific case 
of the Tolman type— IV solution 2], and later by Finch 
and Skea Anisotropic stellar models, with the re- 

spective astrophysical applications, were also extensively 
analyzed in Refs. 0, by considering a specific case of 
the Matese- Whitman mass function. The numerical re- 
sults outlined show that the basic physical parameters, 
such as the mass and radius, of the model can describe 
realistic astrophysical objects such as neutron stars [i^ . 

Using Eq. ®, g(r) is given by 



. X ,'bor 



(1 + 3t^) + (1 + cj)26or-2 



(15) 



(l + 6or2)(l + 25or2) 

which is plotted in Fig. [3 Note that g{r) > 0, for 
Lu > —(1 -I- 2bor'^)/{3 + 25or2), indicating an inwardly 
gravitational attraction; from the plot one verifies that, 
qualitatively, g(r) is positive for values of uj in the neigh- 
borhood of —1/3. Now, to be a solution of a gravastar, 
it is necessary that the local acceleration due to gravity 
of the interior solution be repulsive, so that the region 
above-mentioned for which g(r) > is necessarily ex- 
cluded. 

The spacetime metric for this solution is provided by 

,{l-Lo)/2 



ds^ 



- (1 + bor^ 
1 + 260^2 



dr^ 



[1 + 2bor^ 
■ {d9^ + si 
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1 + 5or2 

The stress-energy tensor components are given by 
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The anisotropy factor takes the following form 

^{(1 + cc;)(3 + 26or') [(1 + -iuj) + 26or'(l + uj)\ 

-8u{5 + 26or2)(l + bor^)}/ [(1 + 2bor^)Hl + bor^)] ,(18) 
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FIG. 2: Plot of the "gravity profile" for the Tolman-Matese- 
Whitman mass function. We have defined the following di- 
mensionless parameters: 7 — g{r)/\/bo and /3 = Vbor. Note 
that, qualitatively, g{r) takes positive values of u in the neigh- 
borhood of —1/3, indicating an inwardly gravitational attrac- 
tion. To be a solution of a gravastar, it is necessary that the 
local acceleration due to gravity of the interior solution be 
repulsive, so that the region for which g{r) > is necessarily 
excluded. See the text for details. 

which is plotted in Fig. O where one verifies A > 0. 
For the particular case of lu — —1, the anisotropy factor 
reduces to 

\Att y (1 + 26or2)3 ■ ^ ' 



tion of state, p = up with oj < —1/3, to an exte- 
rior Schwarzschild vacuum solution with p — p — 0, 
at a junction interface S, with junction radius a. The 
Schwarzschild metric is given by 

ds^ = -(^l-^^ de + + r2 {de^ + sin2 9dq^^), 

(20) 

which possesses an event horizon at — 2M. To avoid 
the latter, the junction radius lies outside 2M, i.e., a > 
2M. We shall show below that M, in this context, may 
be interpreted as the total mass of the dark energy star. 
(In an analogous manner traversable wormhole models 
were constructed in Refs. and the stability analysis 
of phantom wormholes was carried out in Ref. 39] .) One 
may also impose boundary and regularity conditions, and 
their respective constraints 0, ^ , at the center and at 
the surface of the dark energy star, without the presence 
of thin shells, however, this analysis will be presented in 
a subsequent paper. 

Using the Darmois-Israel formalism jZ^, the surface 
stress-energy tensor S** at the junction surface S is given 

by the Lanczos equations, Sj — — g^(Kj- — Sjk'^). is 
defined by the discontinuity of the extrinsic curvatures, 
Kj, across the junction interface, i.e., Kij — — K^. 
The extrinsic curvature is defined as Kij — e^j^eQ-j, 
where is the unit normal 4— vector to S, and e^^.^ are 
the components of the holonomic basis vectors tangent 
to S. Using the metrics Q and H2()(l . the non-trivial 
components of the extrinsic curvature are given by (see 
Ref. js^ for further details) 




FIG. 3: Plot of the anisotropy factor for the Tolman-Matese- 
Whitman mass function. We have defined the following di- 
mensionless parameters: S — A/bo and P = \/6o'"- Note that 
A > for generic to and r. 



and 
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T + 



M 



2M 



m+ujam 
TTi 



(a-2m) 



-, 2m(a) . . 9 
1 — + 



2M 



a V a 



_ 1 / 2m(a) 
a V a 



(21) 



(22) 



(23) 
(24) 



IV. STABILITY OF DARK ENERGY STARS 

A. Junction conditions 

In this article, we shall model dark energy stars by 
matching an interior solution, governed by an equa- 



where the overdot denotes a derivative with respect to the 
proper time, r, and the prime henceforth shall denote a 
derivative with respect to a. Equation evaluated at 
a has been used to eliminate g{r) from Eq. 1)22(1 . In the 
analysis that follows, we shall apply the same procedure, 
so that in deducing the master equation, which dictates 
the stability regions, one shall only make use of the mass 
function 'm{r). Thus, the Lanczos equations provide the 
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surface stresses given by 



I 2M , / 2m , \ , , 



T 




a l~2m/a 



(26) 



a and V are the surface energy density and the tangential 
surface pressure, respectively. 

We shall also use the conservation identity given by 

, where [X\^ denotes the disconti- 



T p 



nuity across the surface interface, i.e., [X]_ — X+js — 
X^Ij]. The momentum flux term in the right hand side 
corresponds to the net discontinuity in the momentum 
flux — Tfj^i, If^ which impinges on the shell. The con- 
servation identity is a statement that all energy and mo- 
mentum that plunges into the thin shell, gets caught by 
the latter and converts into conserved energy and mo- 
mentum of the surface stresses of the junction. 

Note that Sl.^^ = - [& + 2a{a + V)/a], and the flux 
term is given by (see Ref. (^3 for details) 



+ _ {p + Pr)ay^l — 2m/a + i 
- 1 — 2m/a 



(27) 



where p and Pr may be deduced from Eqs. ©-Q, re- 
spectively, evaluated at the junction radius, a. Therefore, 
using these relationships, the conservation identity pro- 
vides us with 



cr' = {<J + V) + E, 



(28) 



where S, defined for notational convenience, is given by 



1 m'(l+uj) / 2m 

7 T A / 1 ^ ^ ■ 

'iira [a — 2m) V a 



(29) 



using Eq. © evaluated at a. Note that the flux term is 
zero when uj = —1, which reduces to the analysis consid- 
ered in Ref. [23. 

Taking into account Eqs. |2Sl)-(ESIl and Eq. then 
Eq. finally takes the form 



1 (l-^+a'-ad 



1 _ 2M ^ ^2 

a 



1 _ 2m _^ ^2 



(30) 



which evaluated at a static solution, ao, shall play a fun- 
damental role in determining the stability regions. (A 



linearized stability analysis of spherically symmetric thin 
shells and of thin-shell wormholes where the flux term is 
zero was studied in Refs. |50l|'). 

The surface mass of the thin shell is given m^ — Aiia^a. 
By rearranging Eq. H25I) . evaluated at a static solution 
ao, one obtains the total mass of the dark energy star, 
given by 



M — m(ao) -f ms(ao) 



'1 - 



2m(ao) ms(ao) 



ao 



2ao 



(31) 



Using TTis = Atto^u, and taking into account the radial 
derivative of cr', Eq. I|28|l can be rearranged to provide 
the following relationship 



Ana'ri , 



(32) 



with the parameter 77 defined as = V' /a' , and T given 
by 



T 



Att , 

— {a + V) 
a 



2'Ka ^ . 



(33) 



Equation (|32|l will play a fundamental role in determin- 
ing the stability regions of the respective solutions, rj 
is used as a parametrization of the stable equilibrium, 
so that there is no need to specify a surface equation of 
state. The parameter ^Jr] is normally interpreted as the 
speed of sound, so that one would expect that < 77 < 1, 
based on the requirement that the speed of sound should 
not exceed the speed of light. We refer the reader to Refs. 
|5nl | for discussions on the respective physical interpreta- 
tion of rj in the presence of exotic matter. 

It is also of interest to analyze the evolution identity, 
given by: [T^^ ^ T ^ where T ^ = {K'+ + 

K'^~)/2. The evolution identity provides the following 
relationship 



Pr + 



{p + Pr)a^ 
1 — 2m/a 




2M 



1 — 2m/a 



1 _ 2m 4, ^2 



(34) 



It is of particular interest to obtain an equation governing 
the behavior of the radial pressure in terms of the surface 
stresses at the junction boundary, at the static solution 
ao, with a — a = Q. From Eq. (|34|l . we have the following 
pressure balance equation 



1 2M 2m 

Pr[ao) = Wl h Wl 

ao \ V ao V ao 




CT.(35) 
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Equation l|35|) relates the interior radial pressure imping- 
ing on the shell in terms of a combination of the surface 
stresses, a and V, given by eqs. if^ - lf^ evaluated at 
the static solution, and the geometrical quantities. To 
gain some insight into the analysis, consider a zero sur- 
face energy density, cr = 0. Thus, Eq. (pl5|l reduces to 



2 /, 2M _ 

Pr(ao) = W 1 V . 

ao V ao 



(36) 



As the pressure acting on the shell from the interior is 
negative Prio-Q) < 0, i.e., a radial tension, then a positive 
tangential surface pressure, 7^ > 0, is needed to hold the 
thin shell against collapse. 



B. Derivation of the master equation 

Equation H25|l may be rearranged to provide the thin 
shell's equation of motion, i.e., + V{a) = 0, with the 
potential given by 



V{a) = Fia) 



ms{a) 


2 


'aG{ay 


2a 




_TOs(a)_ 



(37) 



where, for notational convenience, the factors F{a) and 
G(a) are defined as 

m(a) +M , N M - m{a) 

F(a) = 1 — and G{a) = .(38) 

a a 

Linearizing around a stable solution situated at cq, we 
consider a Taylor expansion of V{a) around ao to second 
order, given by 

V{a) = V{aQ) + V'{ao){a-ao) 

+ ^V"iao){a-aof + O[{a-aof] . (39) 

Evaluated at the static solution, at a = ag, we verify 
that V{ao) = and V^'(ao) = 0. From the condition 
y (ao) = 0, one extracts the following useful equilibrium 
relationship 



r ; 



rus 
2ao 



,(40) 



which will be used in determining the master equation, 
responsible for dictating the stable equilibrium configu- 
rations. 

The solution is stable if and only if V{a) has a local 
minimum at ao and y'(ao) > is verified. Thus, from 
the latter stability condition, one may deduce the master 
equation, given by 



da 



> e. 



(41) 



by using Eq. H32() . where rjo = r]{ao) and 8 is defined by 



crT 



27rao 



(42) 



with 



T 



oG 



-(^)(^^) • (43) 



aG 
m 



Now, from the master equation we find that the stable 
equilibrium regions are dictated by the following inequal- 
ities 



Vo < ^, 
with the definition 

n = Q 



if 
if 



da^ 
da 



da^ 
da 
da^ 
da 



>0, 



<0, 



(44) 
(45) 

(46) 



We shall now model the dark energy stars by choosing 
specific mass functions, and consequently determine the 
stability regions dictated by the inequalities (|44l) - (|45|l . 
In the specific cases that follow, the explicit form of il is 
extremely messy, so that as in 38, 39], we find it more 
instructive to show the stability regions graphically. 



C. Stability regions 

There is only the need to specify the mass function 
m{r), as the "gravity profile" g{r) has been eliminated 
from the stability analysis in this section, by using Eq. 
© evaluated at a. In the examples that follow, we shall 
adopt a conservative point of view, by interpreting that 
is the speed of sound, and taking into account the 
requirement that the latter should not exceed the speed 
of light, i.e., < ?/ < 1, on the surface layer. We shall 
also impose a positive surface energy density, cr > 0, 
which implies m(a) < M . As mentioned above, we shall 
not show the specific form of the functions 2, and O, 
leading to the master equation, as they are extremely 
lengthy. However, the stability regions will be shown 
graphically. 



1. Constant energy density 

Consider the specific case of a constant energy density, 
with the mass function and the gravity profile, which we 
shall include for self-completeness, given by 



m[r) = Ar^ 

Ar{l + w) 



l{r) 



1 - 2Ar^ 



(47) 
(48) 
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with A = inpo/S. We shall impose a positive surface 
energy density, tr > 0, so that m{a) < M. Note that this 
latter condition, with a > 2M, places an upper bound on 
the constant energy density of the star's interior, namely, 
po < 3/(327rM2). 

For the case of m(a) < M, one may prove that 
da'^/da\ao < 0, so that the stability regions are dictated 
by inequality (|45|l . The latter is shown graphically in 
Fig. 01 for the specific case of m{a) — M/2. 





-0.0001 H 



-0.0002 - 
H 



-0.0003 
-0.0004 j 

-0.0005 




ao/M 



FIG. 4: We have defined the dimensionless parameter H = 
M'^da^ /da\ag- Note that for the case of m{a) < M, we have 
H < 0. The latter is shown for the specific case of m(a) — 
M/2. 

Considering the cases of w = — 0.5 and uj = —1.25, the 
stability regions are given by the plots depicted below 
the surfaces in Fig. Note that the stability regions are 
sufficiently close to the event horizon, which is extremely 
promising. For this case, the stability regions decrease for 
decreasing uj, i.e., as the dark energy parameter drops 
into the phantom regime. Note that, qualitatively, as 
m M , the only stability regions that exist are in the 
neighborhood of where the event horizon is expected to 
form. 

The above analysis shows that stable configurations 
of the surface layer, located sufficiently near to where 
the event horizon is expected to form, do indeed exist. 
Therefore, considering these models, one may conclude 
that the exterior geometry of a dark energy star would be 
practically indistinguishable from an astrophysical black 
hole. 



2. Tolman-Matese- Whitman mass function 

Consider the Tolman-Matese- Whitman mass function, 
which represents a monotonic decreasing energy density 
in the star interior, and the respective "gravity profile" , 
given by 



m{r) 



bor^ 



9{r) = ( — 



2(1 + 26or2) ' 

{l + Zuj) + {\ + uj)2bor^ 



(l + 6or2)(l + 25or2) 



(49) 
(50) 




FIG. 5: Plots of the stability regions for a dark energy star 
with constant energy density. We have considered uj = —0.5 
and UJ = —1.25, in the first and second plots, respectively. The 
stability regions are given below the surfaces. Qualitatively, 
one verifies that the stability regions decrease for a decreasing 
dark energy parameter. 



Note that the constant ho may be expressed in terms of 
the mass function, as 6o = 2m[a^(l — 4m/a)]~^. Now, as 
6o is considered to be non-negative, by construction, then 
an additional restriction needs to be imposed, namely, 
a > 4m(a). However, as we are primarily interested in 
the behavior of the model's surface layer where the event 
horizon is expected to form, i.e., a > 2M, we shall only 
analyze the domain < m/M < 1/2. 

For the case of m{a) < M, as in the previous example, 
one may also prove that da"^ /da\aQ < 0, so that the sta- 
bility regions are also dictated by inequality (|45|) . The 
stability regions are depicted in Fig. [HI for the specific 
case of w = —0.5. It is possible to show that the stabil- 
ity regions, for this specific case, are very insensitive to 
variations in to. One may prove that the stability con- 
figurations slightly increase for high m/M, and slightly 
decrease for relatively low values of m/M , for decreasing 
values of the dark energy parameter uj. 

The message that one may extract, as in the previ- 
ous example, is that stable dark energy stars exist with 
a transition layer placed sufficiently close to where the 
event horizon is expected to form, so that the exterior 
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geometry of these stars would be difficult to distinguish 
from an astrophysical black hole. 




FIG. 6: Plot of the stability region for a dark energy star with 
the Tolman-Matese- Whitman mass function, where we have 
considered the case of — —0.5. The stability region is given 
below the surface. See the text for details. 



V. SUMMARY AND CONCLUSION 

Although evidence for the existence of black holes is 
very convincing, a certain amount of scepticism regard- 
ing the physical reality of event horizons is still encoun- 
tered, and it has been argued that despite the fact that 
observational data do indeed provide strong arguments in 
favor of event horizons, they cannot fundamentally prove 
their existence In part, due to this scepticism, a 

new picture for an alternative final state of gravitational 
collapse has emerged, where an interior compact object 
is matched to an exterior Schwarzschild vacuum space- 
time, at or near where the event horizon is expected to 
form. Therefore, these alternative models do not possess 
a singularity at the origin and have no event horizon, 
as its rigid surface is located at a radius slightly greater 
than the Schwarzschild radius. In particular, the gravas- 
tar picture, proposed by Mazur and Mottola ii2L] . has an 
effective phase transition at/near where the event hori- 
zon is expected to form, and the interior is replaced by 
a de Sitter condensate. The latter is then matched to 
a thick layer, with an equation of state given hy p = p, 
which is in turn matched to an exterior Schwarzschild so- 
lution. A simplified model was then proposed by Visser 
and Wiltshire , where the matching occurred at a thin 
shell. It has also been argued that there is no way of dis- 
tinguishing a Schwarzschild black hole from a gravastar 
from observational data pol |. In this work, a generaliza- 
tion of the gravastar picture was explored, by considering 
a matching of an interior solution governed by the dark 
energy equation of state, uj = p/p < —1/3, to an exterior 
Schwarzschild vacuum solution at a junction interface. 
We emphasize that the motivation for implementing this 
generalization arises from the fact that recent observa- 



tions have confirmed an accelerated expansion of the Uni- 
verse, for which dark energy, a cosmic fluid dictated by 
an equation of state given by = p/p with lu < —1/3, is 
a possible candidate. We have analyzed several relativis- 
tic dark energy stellar configurations by imposing specific 
choices for the mass function. The first case considered 
was that of a constant energy density, and the second 
choice, that of a monotonic decreasing energy density in 
the star's interior. We then further explored the dynam- 
ical stability of the transition layer of these dark energy 
stars to linearized spherically symmetric radial pertur- 
bations about static equilibrium solutions. It was found 
that large stability regions do exist, which are located 
sufficiently close to where the event horizon is expected 
to form, so that it would be difficult to distinguish the 
exterior geometry of the dark energy stars, analyzed in 
this work, from an astrophysical black hole. 

The possibility of the existence of dark energy, respon- 
sible for the present accelerated expansion of the Uni- 
verse, has opened up new possibilities in theoretical re- 
search. In this context, by extending the notion of a spa- 
tially homogeneous dark energy fluid, to inhomogeneous 
spherically symmetric spacetimes, such as dark energy 
stars, we point out that the latter have some interesting 
physical properties. Note that recent fits to supernovae, 
CMB and weak gravitational lensing data favor an evolv- 
ing equation of state, with dark energy parameter cross- 
ing the phantom divide ui = —1 |35l l36| . Thus, in a 
rather speculative scenario, one may consider the exis- 
tence of a dark energy star, with an evolving parameter 
starting out in the range — 1 < w < —1/3, and cross- 
ing the phantom divide, a; = — 1. A possible approach 
would be to consider a mixture of interacting non-ideal 
fluids, as it seems that in a cosmological setting the tran- 
sition into the phantom reg ime is physically implausible, 
for a single scalar held [33 • One could also apply vari- 
ations of the approach outlined in Ref. 37], where an 
inhomogeneous spherically symmetric dark energy scalar 
fleld was considered, although the latter model does not 
allow a phantom crossing. Once in the phantom regime, 
the null energy condition is violated, which physically 
implies that the negative radial pressure exceeds the en- 
ergy density. Therefore, an enormous pressure in the 
center may, in principle, imply a topology change, con- 
sequently opening up a tunnel, and converting the dark 
energy star into a wormhole p3ll51| . (However, it is still 
uncertain whether topology changes will be permitted by 
an eventual theory of quantum gravity). It has recently 
been shown that traversable wormholes may, in princi- 
ple, be supported by phantom energy (32i . i33ij . which 
apart from being used as interstellar shortcuts, may in- 
duce closed timelike curves with the associated causality 
violations ,51^ ,52^ . It should be interesting to construct 
a mathematical model illustrating this conversion, i.e., 
dark energy star into a wormhole, by considering a time- 
dependent dark energy parameter, which we leave for a 
future work. Perhaps not so appealing, one could denote 
these exotic geometries consisting of dark energy stars (in 
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the phantom regime) and phantom wormholes as phan- 
tom stars. As emphasized in the Conclusion of Ref. |25j| , 
we would also like to state our agnostic position relatively 
to the existence of dark energy stars, however, it is im- 
portant to understand their general properties to further 
understand the observational data of astrophysical black 
holes. 
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